Nuclear stopping from 0.09^4 GeV to 1.93^4 GeV and its correlation to flow 
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We present a complete systematics (excitation functions and system-size dependences) of global 
stopping and sideflow for heavy ion reactions in the energy range between 0.09^1 GeV and 1.93^4 
GeV. For the heaviest system, Au+Au, we observe a plateau of maximal stopping extending from 
about 0.2A to 0.8^4 GeV with a fast drop on both sides. The degree of stopping, which is shown 
to remain significantly below the expectations of a full stopping scenario, is found to be highly 
correlated to the amount of sideflow. 

PACS numbers: 25.75.Ld;25.70.Pq 



In the last two decades collisions between accelerated 
heavy ions have become a major laboratory tool in large 
scale efforts to learn more about the properties of hot 
and compressed nuclear matter. Nuclei are small ob- 
jects however: from the Bethe-Weizsacker formula de- 
scribing empirical nuclear masses it follows that surface 
and Coulomb effects, that are not present in the interior 
of macroscopic nuclear objects, such as neutron stars, re- 
duce the bulk binding energy per nucleon of 16 MeV to 
about half the value, i.e. 8 MeV. For this reason the 
one-fluid hydrodynamics approaches 0, 0, 0] that were 
favoured in the early days to link experimental observa- 
tions and fundamental properties of nuclear matter, such 
as the relation between pressure, density and tempera- 
ture (Equation of State, EoS) were not trivially justi- 
fied. The first convincing evidence of collective (side) 
flow 0, 0| raised hopes that this connection could be es- 
tablished successfully 0. Meanwhile, flow phenomena 
in heavy ion reactions have been studied extensively Q. 
With the advent in the mid eighties of event simulation 
codes implementing microscopic transport theory || it 
soon became clear that non-equilibrium effects in such 
reactions could not be ignored. In one of the most recent 
attempts to derive constraints on the EoS from heavy 



ion data (which also mentions some of the earlier works) 
it is made clear that progress on this topic requires im- 
proved understanding of the momentum dependence of 
the mean fields generated in a heavy ion collision (first 
introduced in refs. as well as an extensive ad- 

justment to experimental information on the degree of 
stopping achieved. 

This Letter presents a complete systematics of stop- 
ping properties in heavy ion collisions at SIS energies. 
Our data show that the degree of stopping achieved in 
a heavy ion collision 1) reaches a well defined plateau of 
maximal stopping centered around 0.5A±0.3A GeV with 
a fast drop on both edges, 2) stays significantly below the 
hydrodynamic limit at all energies, 3) shows a strong cor- 
relation to the global sideflow suggesting that parametric 
adjustments to both observables might require an itera- 
tive procedure. 

The experiments were performed at the heavy ion ac- 
celerator SIS of GSI/Darmstadt using the large accep- 
tance FOPI detector • A total of 25 system-energies 
are analysed for this work (energies in A GeV are given 
in parentheses): 40 Ca+ 40 Ca (0.4, 0.6, 0.8, 1.0, 1.5, 1.93), 
58 Ni+ 58 Ni ( _ 09; o 15) o 25) o.4), 96 Ru+ 96 Ru (0.4, 1.5), 

129 Xe+ 133 Cs 12 7l ( 15j 25; o 4); 197 Au+ 197 Au (q.09, 
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0.12, 0.15, 0.25, 0.4, 0.6, 0.8, 1.0, 1.2, 1.5). Particle 
tracking and energy loss determinations are done using 
two drift chambers, the CDC (covering polar angles be- 
tween 35° and 135°) and the Helitron (9° - 26°), both 
located inside a superconducting solenoid operated at a 
magnetic field of 0.6T. A set of scintillator arrays, Plas- 
tic Wall (7° - 30°), Zero Degree Detector (1.4° - 7°), 
and Barrel (42° — 120°), allow us to measure the time 
of flight and, below 30°, also the energy loss. For the 
low-energy experiments (up to 0AA GeV) the Helitron 
was replaced by an array of gas ionization chambers fl3j 
allowing charge identification of heavier clusters. The ve- 
locity resolution below 30° was (0.5 — 1.5)%, the momen- 
tum resolution in the CDC was (4— 12)% for momenta of 
0.5 to 2 GeV/c, respectively. Use of CDC and Helitron 
allows the identification of pions, as well as good isotope 
separation for hydrogen and helium clusters in a large 
part of momentum space. Heavier clusters are separated 
by nuclear charge. 

Collision centrality selection was obtained by binning 
distributions of cither the detected charge particle multi- 
plicity, MUL, or the ratio of total transverse to longitudi- 
nal kinetic energies in the center-of-mass (c.o.m) system, 
Erat. 

The scaled directed flow is p^} ir = p x dir/ui cm where 
Pxdir = J2 si 9 n iy) Zu x/J2 z ( z fragment charge, u lcm 
spatial part of the c.o.m. projectile 4-velocity, u x = (3 X ^ 
projection of the fragment 4-velocity on the reaction 
plane ^3)- The sum is over all measured charged 
particles with Z < 10, excluding pions, and y is the 
c.o.m. rapidity. In the left panels of Fig. ^ we show 
the scaled impact parameter, = b/b max , dependence 
of p x °} ir for three system-energies, 40 Ca + 40 Ca at 0.4^4 
GeV, and 197 Au+ 197 Au (OA A and 1.5 A GeV). We take 
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bmax — 1.15(^4p° + ) fm and estimate b from the 
measured differential cross sections for the Erat or the 
multiplicity distributions, using a geometrical sharp-cut 
approximation. The Erat selections show better impact 
parameter resolution for the most central collisions than 
the multiplicity selections. The maximum value of p x °J ir 
is a robust observable that does not depend significantly 
on both the selection method and the apparatus filter. 
For a comparative study of alternate, more differential 
measures of sideflow in the published literature we refer 
to the reviews Q. 

To obtain information on the degree of stopping, we 
measure rapidity distributions of the emitted particles. 
In an XYZ Cartesian coordinate system with the Z-axis 
oriented in the beam direction and the X-axis in the reac- 
tion plane, we define the (conventional) longitudinal ra- 
pidity distribution in the beam direction, and the trans- 
verse rapidity averaged over the X and Y directions by 
projecting on a transverse laboratory fixed axis. All ra- 
pidities arc in the c.o.m. and are scaled by the beam 
rapidity. Results for the same system-energies as in the 
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FIG. 1: Left panels: global sideflow as a function of the 
scaled impact parameter (triangles: Erat selection, crosses: 
multiplicity selection). Right panels: global scaled transverse 
and longitudinal rapidity distributions for the Erat selection 
ft' ' < 0.15. Data for three indicated system-energies are 
shown. Systematic errors are discussed in the text. 



left panels are shown in the right panels of Fig. 

The distributions are full event observables: all emit- 
ted charged baryonic particles are added, weighted with 
their respective nuclear charges. The total charge ac- 
counted for represents typically 95% of the system charge 
or more, except for data below 150 A MeV, where thresh- 
old effects prevent complete coverage of fragments with 
charge Z > 6. Where necessary, bi-dimcnsional fitting 
procedures in u t vs y space, which require reflection sym- 
metry, have been used 15] to interpolate and partially 
extrapolate the measured spectra (u t = [3 t ^f transverse 
4-velocity, y = y z longitudinal rapidity). Extrapolation 
corrections to obtain 47r -distributions were less than 10% 
for protons and clusters with nuclear charge Z = 1, 2 and 
up to a maximum of 30% for heavier fragments. 

As a measure of the degree of stopping we present data 
for an observable dubbed vartl, the ratio of the variances 
of the transverse to that of the longitudinal rapidity dis- 
tributions. For an isotropic thermal source vartl should 
be equal to unity. Flow and transparency may change 



vartl however. Early, ideal hydrodynamics estimates 
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FIG. 2: Excitation functions (left) and size dependences 
(right) of the degree of stopping (upper panels) and maxi- 
mal global sideflow (lower panels). The dashed curves result 
from a correction for statistical fluctuations in the reaction 
plane orientation. The curve labelled 'hydro' is an estimate 
of the hydrodynamic limit for the sideflow. See text. 



predict in this energy regime the phenomenon of squeeze- 
out^ resulting, for central collisions, in enhanced emission 
at polar (c.o.m.) angles of 90°: vartl then becomes signif- 
icantly larger than unity. To avoid the uncertainties and 
distortions of the measured far-out tails of these distri- 
butions, we limit the integration to the interval between 
+1 and -1 in the scaled values. By definition Erat and 
vartl are related, although technically different. Auto- 
correlations in vartl, as well as p x dir, were avoided by 
removing the particle of interest from the selection cri- 
terion Erat or, respectively, the transverse-momentum 
vector determining the reaction plane. 

The observable vartl is also related in spirit to R = 
(2/ 71 ") J2 \vu\l J2 \Pzi\ suggested in ref. and used in 
earlier stopping studies |17j . Longitudinal rapidity distri- 
butions for selected particles, measured with the PLAS- 
TIC BALL [13 were published [3 for a limited set of 
systems and beam energies. 

At this point we note that all statistical errors in this 
work are significantly smaller than the systematic er- 
rors. Estimates of the latter, varying between 5 and 
12 % for the rapidity distributions, are based on compar- 
isons of measured yields in the c.o.m. backward/forward 
hemispheres, checks of the overall charge balance and 
technical variations of the fitting procedure. For p^} ir 
(vartl) a constant uncertainty of 0.008 (0.04) describes 
best the overall error. Point-by-point fluctuations (shown 



in Fig- EJ) are smaller: 0.005 (0.02). Simulations with the 
transport code IQMD (based on Quantum Molecular Dy- 
namics 0), together with detector filters, were also used 
to assess these errors. 

The upper-left panel of Fig. 2 shows excitation func- 
tions of vartl for Au+Au and Ca+Ca in central collisions 
corresponding to scaled impact parameters < 0.15. 
Interesting aspects of these data are a) the occurrence, 
for Au+Au, of a plateau of maximal stopping extending 
from 0.2 A GeV to 0.8 A GeV with a fast drop on both 
sides, b) the fact that vartl never reaches values above 
or even close to unity, and c) the large system-size de- 
pendence evidenced by comparing the Au+Au with the 
Ca+Ca data. 

The smaller values for the lighter system offer a key 
argument in favour of interpreting vartl < 1 as evidence 
for partial transparency, rather than the dominance of 
longitudinal over transversal pressure gradients. 

More detailed information on the system-size depen- 
dence at two incident beam energies (0.4A GeV and 1.5A 
GeV) is shown in the upper right panel of Fig. 2 which 
also demonstrates that, especially at the lower energy, 
there is no evident saturation as one goes from Ca+Ca, 
via Ni+Ni, Ru+Ru, Xe+Csl all the way to Au+Au. 

The two lower panels of the figure show by compari- 
son with the two upper panels an impressive correlation 
with the maximum (see Fig. ^) scaled directed sideflow, 
max {p^dir\> both in the excitation functions and in the 
size dependence. In particular sideflow is positively corre- 
lated with vartl as shown in Fig. 3 for 22 system-energies 
(data for E/A < 0.15GeV were omitted because only 
upper limits on vartl could be determined). 

Although the two global observables are measured at 
different average impact parameters, the most natural 
explanation of the observed correlation is that sideflow 
is resulting from the pressure gradient between partici- 
pant and spectator matter, which in turn grows with the 
degree of stopping. Note that properly scaled sideflow 
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FIG. 3: Correlation between the maximal directed sideflow 
and the degree of stopping. The line is a linear least squares 
fit to the data, which extend from 0.15 to 1.93 A GeV. 
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is essential to reveal this correlation. Both obscrvables 
are also expected to be influenced by t he statistics of ele- 
mentary collisions. The estimated |23| effect on sidcflow 
is shown in Fig. 2. 

Using the Rankine-Hugoniot-Taub one-dimensional 
shock equations together with a realistic equation-of- 
state (SLy230a from ref. (24|), we have estimated the 
maximal densities and pressures V to be expected 
under the full stopping assumption for the energy 
regime of interest here. Characterizing the system size 
by R S ys and the (passing) time by t pass , the scaled 
global sidcflow should behave as gradV x t pass /ui cm ~ 
(V I 'R S y S )(R sys /u lcm ) /uicm = P/wicro- Tnc l attcr quan- 
tity, with a normalization to be discussed below, is plot- 
ted in Fig. El (lower left panel) and compared to the mea- 
sured scaled sideflow. Note that ideal-gas one-fluid hy- 
drodynamics predicts \'2l\ a size independent flat excita- 
tion function (R sys drops out, and V rises linearly with 

u lcm)- 

Our estimate of the shape of the excitation function 
does not yield the absolute value of max [p^ ir ], which 
depends on the geometrical details of the collision. We 
opted for an estimate of the normalization by perform- 
ing an IQMD calculation using, instead of the free-space 
nucleon-nucleon elastic cross sections, values enhanced 
by a factor two over the known free space values, thus 
mocking up approximately the full stopping situation (we 
found vartl = 1.23, confirming qualitatively early predic- 
tions 0). The resulting maximal sideflow value (0.225 at 
0.4^4 GeV) is marked by a full square in the Figure . This 
serves to give an idea on what to expect in the ideal hy- 
drodynamics (non- viscous, one-fluid) limit. Introducing 
still higher elementary cross sections in the simulation 
does not change the prediction of vartl significantly as 
expected in the hydrodynamic limit and confirming the 
relevance of vartl to assess stopping. 

The difference between the theoretical estimate and 
the experimental observation is likely to be the finite 
viscosity of hot nuclear matter, or, equivalently, the fi- 
nite ratio R S y S /Xf (which is the Reynolds number with 
A / the mean free path) , which is not sufficiently large to 
produce fully macroscopic ideal hydrodynamic behaviour 
in such collisions. The significant difference of typically 
30% in the plateau region, together with the observed 
strong system-size dependence (already known from ear- 
lier studies 0,^3 of the flow of selected ejectiles), stress 
the 'small-object' character of nuclei: finite-size (surface) 
effects are large. The most interesting, and new, aspect 
of the full-event data, however, is the rapid drop away 
from the plateau. In particular, the gradual rise of the 
theoretical prediction at the low energy end is totally 
suppressed by both the onset of partially attractive mean 
fields and the increase of Pauli blocking below 0.4A GeV 
[22^ . On the high energy end (1.5A GeV) of our data, 
looking again at Fig. one can estimate that the effec- 
tive average pressures achieved are as much as a factor 



two lower than the full-stop scenario would predict. In- 
medium modifications of elementary cross sections |25j . 
caused possibly by nucleon-mass changes at high density, 
might be needed to account for the apparent rise of trans- 
parency in an energy regime where strong effects of Pauli 
blocking are no longer expected. 

Our observation of partial transparency confirms qual- 
itatively earlier isospin tracer experiments and, to- 
gether with its close correlation to sideflow, has im- 
portant consequences for the ongoing efforts to extract 
the Equation of State from heavy ion data in the non- 
equilibrium situation |2fij . as well as for attempts to un- 
derstand in-medium effects on mcsonic and baryonic par- 
ticles, since the effective densities in the early phase of 
the reactions are expected to be strongly correlated to 
the degree of stopping. 

This work was partly supported by the German BMBF 
under contract 06HD953 and RUM-99/010 and by KRF 
under Grant 2002-0 15-CS0009. 
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